subroutine set_param
  use param
  use global
  implicit none
  real(dbl) :: t1,t2
  ! ---- initialization ---- !
  ! -- mode
!  rmode=0
  ! -- dise
!  dise=-1
  ! -- sube
!  sube=-1
  ! -- xsec
!  xmode=0
!  xmine=-1
!  xmaxe=-1
!  xeh=-1
!  xih=-1
  ! -- power
!  pmode=0
!  pmine=-1
!  pmaxe=-1
!  peh=-1
!  pih=-1
!  nhp=100
!  nip=10
  ! -- subcascade
!  next_e=0
!  int_e=0
  ! ---- Kernel data ---- !
  e2=(2.2d+6*ph)/(2.0*pi)
  mtoev=9.31494028d+8/(c**2.0)
  kl=2.92
  km=0.191
  kq=0.512
  kex=2.0/3.0
  ikq=1.0/kq
  ka=(0.885*a0)/((2.0*acharge*kex)*0.5)
  k1=(ka/(2.0*acharge*acharge*e2))**2
  k2=0.5*pi*ka**2
  k3=(ka/(2.0*acharge*acharge*e2))**2
  t1=8.0*pi*e2*a0*acharge*acharge**(7.0/6.0)
  t2=(amass*mtoev)**0.5*v0*(2*acharge**kex)**1.5
  k4=t1/t2
  ! ----
  ngi=10
  ! ----
  xsec%ndata=0
  xsec%pos=0
  nstop%ndata=0
  nstop%pos=0
end subroutine set_param
